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Using the Color-Dielectric model to describe quark confinement, including strange quarks and 
accounting for beta-equilibrium, we get masses for a static neutron star in the range f .3 < M/Mq < 
f .54 for a radius R sa 9 km. Rapid cooling through the direct URCA mechanism is obtained. The 
implications for the composition of neutron stars is discussed. 

PACS numbers: 97.60. Jd, 12.39.-x, 24.85.+p 

The equation of state (EOS) for dense matter is central to calculations of neutron-star properties, such as the mass 
range, the mass-radius relationship, the crust thickness and the cooling rate, see e.g. Refs. [[l|-| . The same EOS is 
also crucial in calculating the energy released in a supernova explosion. 

The typical density range of a neutron star stretches from central densities of the order of 5 to 10 times the nuclear 
matter saturation density no = 0.17 fm -3 to zero at the edge of the star. Clearly, the relevant degrees of freedom 
will not be the same in the crust, where the density is much smaller than no, and in the center of the star where the 
density is so high that models based solely on interacting nucleons are questionable. Data from neutron stars indicate 
that the EOS should probably be moderately stiff in order to support maximum neutron star masses in a range from 
approximately IAMq to l.SLM© [0]. In addition, simulations of supernovae explosions seem to require an EOS which 
is even softer. A combined analysis of the data coming from binary pulsar systems and from neutron star formation 
scenarios can be found in Ref. |5[ , where it is shown that neutron star masses should fall predominantly in the range 
1.3 < M/Mq < 1.6. 

Several theoretical approaches to calculations of the EOS have been considered. The hypothesis that strange 
quark matter may be the absolute ground state of the strong interaction ||, has been used by several authors in 
the investigation of the possibility of interpreting pulsars as rotating strange stars . Other approaches introduce 
exotic states of nuclear matter, such as kaon P-|l2[ or pion condensation Jl3| , pT| . 

The scope of this work is to study properties of neutron stars like cooling rate, total mass and radius employing a 
massive quark model, the so-called Color-Dielectric model (CDM) [^5|-^7|. A very important feature of matter, as 
described using the CDM model, is that the deconfinement phase transition happens at a density of the order of 2-3 
times no, a density expected to be reached in the interior of a neutron star. Neutron stars, when studied using the 
CDM, consist of an exterior region described, e.g., by the equation of state of Friedman and Pandharipande jl|], and 
an interior region of quark matter. Moreover, the effective quark mass in the CDM is always larger than a value of 
the order of 100 MeV, hence chiral symmetry is broken and the Goldstone bosons are relevant degrees of freedom. 
This is to be contrasted with models like the MIT bag, where quarks have masses of a few MeV. We therefore expect 
the CDM to be relevant for computing the cooling rate of neutron stars via the URCA mechanism, as suggested by 
Iwamoto [fl9f . 

The CDM is a confinement model which has been used with success to study properties of single nucleons, such as 
structure functions |2(J and form factors |pl[ , or to describe the interaction potential between two nucleons p2j, or 
to investigate quark matter Jl7],^3|]. In particular, it is possible, using the same set of parameters, both to describe 
the single nucleon properties and to obtain meaningful results for the deconfinement phase transition Jl~7| | . The aim 
here is therefore to see whether this model gives reasonable predictions for neutron star observables. 

The Lagrangian 1 of the model is given by: 

C = iftfW 

+ V" 4^~^l (f + Z7 5 T • 7?) "0/ + — ■Ips'lps 

+ \(d,xf-V( X ) (1) 



Throughout this paper we set G = c = % = 1, where G is the gravitational constant. 
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where U(a,ir) is the "mexican-hat" potential, as in Ref. pifl . The lagrangian C describes a system of interacting u, 
d and s quarks, pions, sigmas and a scalar-isoscalar chiral singlet field x whose potential U (x) is given by 

V( X ) = \M 2 x 2 . (2) 

The coupling constants are given by g U) a = g(fn ± £3) and g s = <?(2 fx — fir), where / w = 93 MeV and fx = 1 13 
MeV are the pion and the kaon decay constants, respectively, and £3 = f K ± — f K o — —0.75 MeV. These coupling 
constants depend only on a single parameter g. This parameter is fixed, together with the mass M. of the x ne ld so as 
to reproduce experimental nucleon mass and radius. The SU(3)f version of the model has been introduced by Birse 
and McGovern pSpl . 

When considering a single hadron, confinement is obtained via the effective quark masses m u ^ = g Ui dO I '(x/jt) and 
m s = g s /x, which diverge outside the nucleon. Indeed, the classical fields x an d a are solutions of the Euler-Lagrange 
equations and x g° es asymptotically to zero at large distances. 

The model parameters g and M. are fixed to reproduce single nucleon properties. In order to describe the single 
nucleon state, a double projection on linear and angular momentum eigenstates has to be performed, see Ref. |24| . We 
will use the parameters g = 0.023 GeV and M. — 1.7 GeV, giving a nucleon isoscalar radius of 0.80 fm (exp.val.=0.79 
fm) and an average delta-nucleon mass of 1.129 GeV (exp.val. =1.085 GeV). A similar set of parameters has been 
used to compute structure functions |^0| and form factors |^l[ . 

The quark matter (QM) phase is characterized by a constant value of the scalar fields and by using plane waves to 
describe the quarks. The total energy of QM in the mean field approximation reads 

f=u,d,s J 1 1 

+ VU( X )+VW{a,% = 0), (3) 

where k F is the Fermi momentum of quarks with flavour /. 

The high-density matter in the interior of neutron stars is described by requiring the system to be globally neutral 

(2/3K - (l/3)n d - (1/3K - n e = 0, (4) 

where n u ^ s e are the densities of the u, d and s quarks and of the electrons, respectively. Morover, the system must 
be in /3-equilibrium, i.e. the chemical potentials have to satisfy the following equations: 

Aid = Mu+A*e, (5) 

Ms = Mu + A*e- (6) 

Eqs. (§)-(^|) have to be solved self-consistently together with field equations, at a fixed baryon density n = n u +nd+n s . 

In the upper part of Fig. [I], the resulting quark masses are shown as function of density. As can be seen, the quark 
masses are large and of the order of 100 MeV up to densities about 10 times no- In the lower part of Fig. [[J the Fermi 
momenta of the quarks are shown. In the relevant region for neutron stars the momenta are of the order of 400 MeV. 
Strange quarks appear at a density of ~ l/3no. 

The energy of the QM, as described by the CDM, is larger than the energy of nuclear matter for densities smaller 
than 2-3 times uq. For lower densities we take the results from Friedman and Pandharipande |Q, listed in Table 5 
of Ref. @. We obtain an EOS valid at all densities by joining the EOS of flf] and the CDM EOS through a Maxwell 
construction. The latter starts at a density ~ O.lfm -3 and terminates at a density ~ 0.26fm~ 3 . The energy per 
baryon E/A in beta-stable nuclear matter is displayed in the upper part of Fig. || as function of density n. In the 
lower part of Fig ^ we show the above Maxwell construction. For the sake of comparison we also include the results 
from pure neutron matter obtained with the so-called Walecka model [^7j and the energy per particle obtained from a 
relativistic Dirac-Brueckner-Hartree-Fock (DBHF) calculation in neutron matter p8| , using modern meson-exchange 
potential models ]29|] . As can be seen from Fig. ||, our CDM EOS is much softer than those obtained with either the 
Walecka model or the DBHF calculation. 

From the general theory of relativity, the structure of a static neutron star is determined through the Tolman- 
Oppcnheimer-Volkov equation, i.e. 

dP _ {p{r) +P(r)} {M(r) +4nr 3 P(r)} 

~d7 ~ r 2 - 2rM(r) ' ^' 
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where P(r) is the pressure and M(r) is the gravitational mass inside a radius r. The pressure is related to the energy 
per particle £ = E/A by the relation P(n) = n 2 (d£/dn). 

The resulting mass is shown in Fig. |[ while the total mass as function of radius is displayed in Fig. ||. With 
the CDM we obtain a maximum mass M max m 1.54Mq and a radius of 9 km at a central density corresponding to 
approximately 9 times nuclear matter saturation density, in good agreement with the experimental values for the 
mass [^j^] and the estimate for the radius |p0|| . Concerning the minimum mass, we obtain a value equal to 1.3M Q . 
It is also possible to compute the minimum period of rotation of the star, given by P m ~ (O.8ms)i?^Q 2 (M0/M) 1 / 2 , 
where Rio is the radius in units of 10 km In our model we get P m = 0.55 ms. 

Finally, we study the neutrino and antineutrino emission, which can be responsible for rapid cooling of neutron stars 
(URCA mechanism). In the interior of the star, where, in our model, strange quark matter is present, the relevant 
reactions are: 

d — > u + e~ + v e , (8) 
e~ + u — » d + v e , (9) 
s — » u + e~ + v e , (10) 
e~ + u — » s + V e . (11) 
To conserve momentum in the reactions, the following inequalities have to be satisfied: 

\kp kp\ < kp kp -\- kp, (12) 
\kp kp\ kp < kp -\- kp. (13) 
For densities larger than 0.53 fm~ 3 , conditions (112) are satisfied and the direct URCA mechanism involving only u 



and d quarks can start. Conditions (13) are satisfied only for densities larger than 1.4 fm -3 . The URCA mechanism 
involving strange quarks is suppressed by a factor sin 2 9 C ~ 0.05 with respect to the previous process, where 9 C is the 
Cabibbo angle. 

We compute the neutrino and antineutrino luminosity, using the ultrarelativistic expansion for the chemical poten- 
tials ]l9[| . This approximation is reasonable because the ratio m q /kp is approximately equal to 0.25. One gets for the 
total luminosity |19t | 

457 

-Gpcos 2 e c m 2 d fk F (k B T) 6 , (14) 
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where / = 1 — (m u /md) 2 {kp/k F ) — (m e / md) 2 (k F / kp) . To obtain the characteristic cooling time we equate the energy 
loss per unit volume to the rate of change of thermal energy per unit volume r = cyT/e. Here the heat capacity of the 
QM is cy = J2f= u d s mjkpkgT. We get r = C lday/Tg , where T g is the temperature measured in units of 10 9 °K 
and C is a constant ranging from 0.5 to ~ 10 going from heavy to light neutron stars. To compare with experiment, 
one should also account for possible reheating mechanisms. There are however indications ]32] | that temperatures of 
young (~ 10 4 years old) neutron stars lie below that obtained through the indirect URCA processes. This suggests 
that cooling through exotic mechanisms, as accounted for in our model, may be a meaningful explanation. Moreover, 
a recent reappraisal of neutrino-pair bremsstrahlung by Pethick and Thorsson J|3| indicates that this process is much 
less important for the thermal evolution of neutrons stars than suggested by earlier calculations. In addition, since the 
results of Ref. |l|,|^] indicate that the crust of the neutron star is considerably smaller than previously extimated, and 
the results of Ref. [Q show that proton superconductivity in the stellar interior is substantial, thereby inhibiting the 
traditional URCA mechanisms, there is a possibility that exotic states of matter like QM provide a viable explanation 
for various properties of neutron stars. 

In summary, here we show that the CDM, which has been successful in describing various properties of hadrons, also 
gives reasonable results for neutron star observables like the mass and radius. Further, our results for the cooling of 
the star, combined with the results of Refs. may suggest that the core of the neutron star can be described 

by a QM phase. Lastly, supernovae simulations require an EOS which is soft like ours. 
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FIG. 1. The upper part shows quark masses as function of density n. In the lower figure we give the quark composition as 
function of density n. In both figures, solid lines are the results for u quarks, dashed lines correspond to d quarks and dotted 
lines are the results for s quarks. 
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FIG. 2. The upper part gives the energy per baryon E/A as function of the density n for the three equations of state 
discussed in the text. The solid line corresponds to the result obtained with the CDM, the dotted line represents the result for 
the Walecka model, while the dashed line is the result from the DBHF calculations of Ref. [28]. In the lower figure, we show 
the Maxwell construction discussed in the text. The solid line is the final EOS after the Maxwell construction, the dash-dotted 
lines to the left (low densities) is the EOS of Ref. [18], while the dash-dotted line to the right is the CDM EOS. 
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FIG. 3. M/Mq as function of n/no, where no is the nuclear matter saturation density no = 0.17 mi -3 . The solid line 
corresponds to results obtained with the CDM, the dotted line represents results for the Walecka model, while the dashed line 
are results from the DBHF calculations of Ref. [281. 
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FIG. 4. Total mass M/Mq as function of R. Notations are the same as in Fig. 3. The arrows indicate the maximum and 
minimum mass obtained with the CDM EOS. 
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